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Abstract. 

The recombination model as a model for hadronization from a quark-gluon plasma 
I has been recently revived since it has advantages in explaining several important 

features of the final state produced in heavy- ion collisions at RHIC, such as the 
constituent quark number scaling of the elliptic coefficient versus the transverse energy 
of identified hadrons, the bending shape of the px spectrum of hadrons near 5 GeV/c, 
and the measured large value of baryon to meson ratio (of the order of unity) in 
O ■ the same pt range. We have developed a dynamic simulation model of heavy-ion 

' collisions in which a quark-gluon plasma, starting from a certain initial condition, 

evolves hydrodynamically until it reaches the phase boundary, and then hadronizes 
by valence quark recombination. Rescattering after hadronization is described by 
^ ' UrQMD. We discuss some details of the model and report first, preliminary results. 



1. Introduction 



^ ■ The quark recombination model of hadron formation, which was originally proposed in 
1970s[l] and studied actively in 1980s[2j, has been revived|3l HI El [3, [8] as it naturally 



explains some of the features of the produced hadrons in relativistic heavy-ion collisions. 
The hydrodynamic analysis of the measured elliptic flow coefficient of hadrons, which 
signals the collective behavior of the produced system, could explain the magnitude 
and the ordering of the magnitude of the coefficient V2 of various hadrons and thus 
showed that the system produced during the heavy-ion collisions behaves as a strongly 
interacting, nearly perfect fluid. If one plot the elliptic coefficient per constituent quarks, 
V2/n versus the transverse kinetic energy per constituent quarks, Et/ti, then the data 
for various hadrons converge into a single line, suggesting that the collectivity is the 
property not of the hadrons but also of the constituent (valence) quarks. A natural 
explanation is that the hadrons are produced by recombining the quarks which are 
collectively moving in the quark-gluon plasma. 

A second set of data supporting the recombination of quarks to form hadrons 
are the transverse momentum spectra of hadrons, which shows bending behavior near 
~ 5 GeV. Hadrons with small momentum [pt ~ 2 GeV) are thermal ones and 
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large momentum hadrons are those from jets. Fragmentation of jets may produce 
hadrons with pr ^ 5 GeV while the recombination of quarks with thermal distribution 
in the quark-gluon plasma can produce the hadrons with pr ^ 5. The two mechanisms 
combined have shown to fit the hadron spectra in the mid-p^ range well. Another 
experimental observation supporting the recombination model is the large baryon-to- 
meson ratio, e.g. R{p/tt) ~ 1, near pj' ~ 2 GeV. If the hadrons are formed by 
recombining quarks with a thermal distribution, the probabilities of forming mesons 
or baryons are determined only by the probabilities to find quarks, and the ratio 
R{p/tt) ~ 1 can be naturally explained. 

The recombination calculations reported so far were mostly done in the framework 
of schematic evolution models. The only dynamical calculation is that by Ko et. 
ai: The initial partonic matter evolves through the parton cascade and at the phase 
boundary hadrons are formed by recombination. Thus fixing the initial state just after 
the the nuclear impact, they can simulate the main experimental features described 
above. Here we report a dynamic simulation method describing the heavy-ion collisions 
where the evolution of the quark-gluon plasma is described by ideal hydrodynamics 
up to the phase boundary. At the phase boundary, hadronization is simulated by the 
recombination mechanism and the produced hadrons may or may not further rescatter 
until thermal freeze-out. In our implementation of this model, quark-gluon evolution 
is described by the hydrodynamic equations. For the hydrodynamic cells at the phase 
boundary, hadronization occurs via recombination model and thus the Monte Carlo 
program creates hadrons with their definite position and momentum at each cell whose 
expansion velocity is explicitly known. The hadrons may further undergo hadronic 
rescattering and this process is described by UrQMD. 

2. Dynamic Recombination Model 

We use the three-dimensional ideal relativistic hydrodynamic program by C. Nonaka [9] 
which uses comoving (Lagrangian) grid points rather than the commonly used fixed(in 
some coordinate system) Eulerian grid. The advantage of the Lagrangian grid for the 
calculation the hypersurface integral of the type appearing in the Cooper-Frye formula 
will be explained below. 

Each quark-gluon plasma fluid cell at Tc is specified by the temperature, Tc, 
the baryon chemical potential, fiB, position, {T,x,y,ri), and the spatial length, 
{dr, dx, dy, dr]), where r = \Jt'^ — is the proper time and ri = artanh(2;/t) is the 
space-time rapidity. The cell is moving with its 4-velocity = 7(1,-1?). The quark 
distribution in the cell at hadronization is given as 




(1) 



The invariant cross sections for recombined hadrons are given as [5]: 




(2) 
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gM(P) ^ j^P^da^j d^k \(l)M(x)\^Wa(k)wb(P - k) (3) 
for mesons, and 

„ dNn - Cj , „s , , , 

9b{P)^ I P^'da^, j d^kid^k2\(t>B{xuX2)\^Wa{ki)Wi{k2)Wc{P-ki-k2){f>) 

for baryons, where S denotes the integration over the hadronization hypersurface. 
Inserting the quark distribution into the above equations, one gets 

gM{P)^ j p^da^eMP-^/T) J d^keM-l(Ea + E,)/T)\(l>M(k)\'' (6) 

for mesons, where Ea and E^ are the energies of the quark and antiquark, respectively. 
Similarly for baryons: 

gsiP)^ jp^da^eMP-'^/T) j d''kA\Mki.k2)\^ 

xeM-liEa + Eb + E,)/T). (7) 

Since the wave functions for hadrons are best known in the light-front frame, 
it is convenient to use the light-cone variable, x, which denotes the longitudinal 
momentum fraction carried by the quark in the infinite momentum frame. If one 
neglects the transverse momentum of quarks relative to the hadron momentum P, the 
equations become simplified. When the wave functions for hadrons are known, one can 
explicitly decompose the quark momentum into the longitudinal transverse components, 

— * — * — * 

ki = xP + ki^±, and evaluate the integrals exactly. In the above equations, integration 
over the hadronization hypersurface is done numerically. The hypersurface areas 
are given by the cell surfaces between adjacent cells, one of which has a temperature 
T >Tc and the other has T <Tc. In the case of a first-order phase transition with phase 
coexistence, also pairs of cells with a different fraction of quark-gluon plasma contribute. 
A single grid cell may have space-like as well as time-like surfaces, contributing to the 
hadronization hypersurface. It is known that the hydrodynamic programs using an 
Eulerian grid sometimes give negative values of the expression P^da^. However, in the 
hydrodynamic program which uses the Lagrangian coordinate such as developed by C. 
Nonaka and is used in this calculation, P^da^ is positive definite. 

During the hadronization energy and momentum should be conserved, equivalent 
to energy conservation in the rest frame of a cell. We enforce energy conservation by 
continuing the creating process of hadrons, using the Monte Carlo method, until all the 
energy of the cell is exhausted for hadronization. However, since our grid is so fine that 
the total energy of each individual fluid cell is smaller than the energy required to create 
a hadron, we add the energies of all fluid cells within a small interval and randomly 
create hadrons from this ensemble of cells. The Monte Carlo implementation of this 
scheme works as follows. Consider a fluid cell at T^. in the mixed phase whose degree of 
hadronization is characterized by the QGP volume fraction a = Vq/{Vq + Vh)- If Ei 
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is the total energy in the fluid cell i and Sar is the negative of the change of a during 
the evolution in one time step, then Ei 6ar is the energy used for the hadronization 
during the time interval. Summing the energy for all the cells in a given 77 interval, one 
flnds the total energy to hadronize among the fluid cells in (77,77 + dr]). Hadronization 
continues until all the energy has been used to create hadrons. 

In the sudden recombination model the probability to create a hadron with certain 
momentum P depends only on the number and flavor of its constituent quarks but is 
independent of the hadron mass. In order to implement these constraints, we dived the 
mesons into four groups with different flavor contents: 

Ml = (qq), M2 = (qs), M3 = (sq), M4 = (ss). (8) 

Similarly, there are a total of eight groups for baryons, four for baryons and four for 
antibaryons: 

Bl = (qqq), B2 = (qqs), B3 = (qss), BA = (sss), 

ABl = (qqq), AB2 = (qqs), AB3 = (qss), ABA = (sss). (9) 

To obtain the probability to hadronize into a hadron with certain flavor content and 
momentum, one randomly choses a cell among those in the same r] interval, the flavor 
group, the momentum and a random number in the interval [0, 1]. When the random 
number exceeds the probability to create the hadron given by the recombination model 
equations, a hadron with certain momentum is created in that fluid cell. This procedure 
is repeated until all the energy available is exhausted. Finally, in order to determine the 
hadron species inside a given group (e.g. p versus A+, or some other resonance included 
in the UrQMD model), we assume a relative Boltzmann weight, die~"^^^'^, where di is 
the degeneracy and rrii the hadron mass. 

As a result of hadronization, we create hadrons moving with a thermal velocity 
spectrum of a fluid cell superimposed on the collective flow velocity of the cell given by 
the hydrodynamics. Rescattering among those hadrons is then simulated with the help 
of UrQMD. By running UrQMD code, we may check whether the chemical freeze-out 
coincides with the time of hadronization, and we can track the cooling of the system 
until the flnal thermal freeze-out. 

3. Results and Summary 

In Fig. [T]the transverse momentum spectra, dN/pfdpt vs. pt of pions and protons are 
shown. Those particles are hadronized from a QGP with an initial state same as in Ref. 
[9] for the U+U collisions at 200 GeV/A, and experienced the flnal state interaction 
by running UrQMD. The results are obtained from 20 events and in order to make the 
systematic error smaller one has to repeat the simulation for many events. In Fig. [2] the 
rapidity spectra of pions are shown. 

As emphasized before, hadrons are generated from a QGP according to the sudden 
recombination model, and rescattering among hadrons is simulated by UrQMD. The 
creation of hadrons with their position and momentum enables us to "measure" any 
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Figure 1. Transverse momentum spectra of tt" (left) and protons (right) hadronized 

from a QGP with the initial state in Rcf.[9], which was used to simulate the central 
Au+Au collisions at 200 GeV-A at RHIC, and experienced the UrQMD afterburner. 



'ubap.11' using 1:4:5 ^ 



Figure 2. tt rapidity spectrum obtained as in Fig. 1 



physical observables of interest. The only inputs are those of the initial state of the 
hydrodynamic evolution and the hadronization temperature Tc (or the equation of state). 
The comparison of the result with the statistical model will be interesting. In the future 
we hope to calculate of the elliptic flow coefficients of various hadrons in order to study 
how the flow anisotropy is related to the initial momentum anisotropy of the QGP. 
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